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Abstract 

Single photon detectors are fundamental tools of investigation in quantum optics and play a central role in 
measurement theory and quantum informatics. Photodetectors based on different technologies exist at optical 
frequencies and much effort is currently being spent on pushing their efficiencies to meet the demands coming 
from the quantum computing and quantum communication proposals. In the microwave regime however, a single 
photon detector has remained elusive although several theoretical proposals have been put forth. In this article, we 
review these recent proposals, especially focusing on non-destructive detectors of propagating microwave photons. 
These detection schemes using superconducting artificial atoms can reach detection efficiencies of 90% with existing 
technologies and are ripe for experimental investigations. 
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1. Introduction 

In 1905, his annus mirabilis , Einstein not only postulated the existence of light quanta (photons) while explaining 
the photoelectric effect but also gave a theory (arguably the first) of a photon detector [1]. In the decades that 
followed, significant progress were made in designing photon detectors based on several technologies such as 
photomultiplier tubes, avalanche photodiodes and cryogenic detectors among others. Such detectors are routinely 
used in experimental setups in applications ranging from spectroscopy to sensors. Recently, there has been a 
huge drive coming from the field of quantum information processing to push the efficiency of photon detectors to 
work at the quantum limit. An ideal photon detector is expected to have 100% efficiency, a very low dark count 
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and a number resolving nature. Such high efficiency single photon detectors are crucial in the implementation 
of quantum cryptography proposals such as quantum key distribution (QKD) [2], in experimental probing of the 
foundations of quantum mechanics such as Bell tests [3] and in implementing all optical quantum computers [4] 
among several other applications. Several single photon detectors have been realized in the optical regime [5,6] 
and are part of the standard quantum optics toolkit. However, single photon detectors at microwave frequencies 
have been difficult to implement and researchers have resorted to special schemes for homodyne and correlation 
measurements [7,8]. The focus of this article is to review some of the proposals that have been put forth to fill 
this gap. 

Of particular interest among photodetectors are the non-destructive ones. These detectors are transparent to the 
incoming photons and the measurement scheme is termed as quantum nondemolition (QND). QND measurements 
were first proposed to detect gravitational waves by evading the measurement back-action on the system [9-13]. 
These measurement schemes were well suited for the field of quantum optics leading to successful implementation 
of QND measurements of photon flux in the optical regime [14]. QND detectors play a major role in schemes such 
as quantum error correction [15], state-preparation by measurement [16,17] and one way quantum computing [18]. 
We will focus on recent proposals for nondestructive detection of microwave photons in section 3. 

We have to nevertheless note that, photon detection schemes (including QND detection) have been shown earlier 
for microwaves stored in cavities [19-22]. While strong interactions could be mediated between the photons in 
the cavity and matter used as a detector, the use of cavities also complicate the setup as one has to worry 
about the bandwidth of operation and the compromise between quality factor and reflection. Schemes to catch 
microwave photons in a cavity by tuning their coupling to transmission line have been developed recently [23-25]. 
However, one has to know the exact shape, width and the arrival time of the wave-packet to completely absorb 
the photon in to the cavity. These additional problems make such schemes less attractive for applications such as 
those discussed above. Hence for the purpose of this review, we will focus on detectors designed for propagating 
microwave photons. 

The current article is organized as follows. In section 2, we will discuss proposals for microwave photon detection 
that are destructive in nature. These proposals are based on current biased Josephson junctions (CBJJ). In the 
next section, we will look at proposals that use the photon-photon interaction mediated using a superconducting 
artificial atom to perform QND detection of microwave photons. In the final section, we will summarize the main 
messages from these proposals. 


2. Photon detectors based on Josephson junctions 

Initial proposals for detecting propagating microwave photons were based on current biased Josephson junctions 
(CBJJ) [26-28]. These junctions have a potential energy U(S) = —(/ c 4>o/27r) cos(J) — /&J, where I c is the critical 
current of the junction, 4>o = h/ 2e is the flux quantum and 5 is the superconducting phase difference across the 
junction. Such a potential, known as a tilted washboard potential (see Fig.1(a)), has several local minima or 
potential wells that have discrete set of states. We can tune the potential barrier of the local minima and the 
number allowed energy levels using the bias current lb . For a particular value of lb (depending on / c ), each of 
the wells acts as a two level system that can be used as a qubit. The lifetime of these levels |0) and |1) depend 
on the tunneling distance to the continuum of modes to the right. As seen from the figure, the tunneling from |0) 
is exponentially suppressed compared to the decay from the state |1). If we neglect the tunneling from the state 
|0), the setup can be mapped to a three level A system (Fig.1(b)), with |1) being a metastable state with a short 
lifetime set by the rate T. The incoming photon with frequency close to the transition frequency uj excites the 
qubit to state |1). The excitation from |1) decays irreversibly with the rate T into the continuum that is labeled 
here as the state | g). This leads to a voltage drop across the junction that can be measured classically and the 
excitation is lost making this a destructive scheme. 

The performance of the detector is characterized by the probability of absorption of the incoming photon. This 
is theoretically calculated as follows. We start with the non-Hermitian Hamiltonian of the A system in the real 
space representation (h = 1) 

H = (uj - iT/2) |1) (1| + iv g J dx (^p^dx^L - ^dx^R^J+V J dx 5(x) [{^ R + ^l) |1) (0| + h.c] , (1) 

where 'ipR/L is the radiation field traveling right/left with a group velocity v g . The interaction between the photons 
and the A system is modeled as a delta potential of strength V at x = 0. A general single excitation wave-function 
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Figure 1. (a) Tilted washboard potential of a current biased Josephson junction. The bias current /{, is chosen such that only two 
energy levels exist within each well, (b) A three level A system that maps to the spectrum shown in (a). 

of the total system is of the form 

\4>) = J dx | ^ R (x,t)^ R (x) + £l(x, |0, vac) + e(t) |l,vac) (2) 

where |0,vac) is the state with the A system in state 0 and the radiation field in vacuum. The single photon 
wave-packets are given by £r/l(x, t). Solving the Schrodinger equation for the above Hamiltonian with this state, 
we get three coupled equations for the coefficients (#,£), ^L(x,t) and e(t). From the solution of these coupled 
equations the probability of level | g) can be calculated as P g m 1 — ||0|| 2 . The value of this probability at long 
times is taken as the measure the detector efficiency. With the A system in the middle of an open transmission 
line, the maximum attainable efficiency was calculated to be 50% [26]. This could be improved by having several 
such scatterers along the line which also increased the bandwidth of operation. In [28], it was theoretically shown 
that by placing one atom at the end of transmission line (in front of a mirror), the efficiency could be increased 
to 100%. 

The above discussed analysis were done for an ideal A system, of which the CBJJ is only an approximation. 
While the tunneling from the state |0) is much smaller than the tunneling from the state 11), it is not zero. 
The direct tunneling from the ground state |0) leads to dark counts, which affects the performance of the setup. 
The performance of CBJJ as a photodetector (referred as Josephson photomultiplier (JPM) in the following 
references to distinguish from phase qubits which ideally work at different parameter regime) was analyzed both 
experimentally and theoretically including the direct tunneling from the state |0) in references [29-32]. However, 
these works are primarily aimed at measuring the photon occupation number inside a cavity. Further analysis is 
needed to check the efficiency of such a system to be used as a detector for itinerant single photons. 


3. QND detection of itinerant microwave photons 

3.1. A transmon in a transmission line 

The transmon [33], a widely used superconducting artificial atom, is a single cooper pair box shunted with a 
capacitance Cs- The large shunt capacitance reduces the charging energy Ec of the circuit that also reduces the 
sensitivity of the qubit to charge noise and is one of the main reasons for the transmon’s popularity. Recently, 
it was shown that a transmon could behave as a Kerr medium that imparts a giant conditional phase shift on 
microwave fields incident on the atom [34]. Following schemes for QND detection of photons based on the cross- 
Kerr effect in the optical regime, proposals for nondestructive single photon detection in the microwave frequencies 
using the transmon were analyzed in references [35-37]. We will briefly review these setups in this section. 

For the purpose of this article, we will consider the transmon as essentially a three level system with the 
Hamiltonian (h = 1) 

H = —ljoi |0) (0| + CJ 12 |2) (2| , (3) 

where ujij — Ej — Ei is the energy difference between levels i and j. The 1 — 2 transition is driven by a coherent 
probe field with frequency u p and amplitude a p . The 0 — 1 transition is driven by the control field of frequency 
uj c which contains either 0 or 1 photon (See Fig. 2). In the rotating frame of the input fields, the Hamiltonian 
can be written including the coherent drive as 

H = -A c |0) (0| + A p |2) (2| + Q p (L 12 + L 21 ), (4) 
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Figure 2. A scheme for single photon detection using a single three level system. Such a setup could be achieved in superconducting 
circuits with a transmon at the end of a transmission line. The difference in the homodyne current with and without a single control 
photon in the wave-packet £(t) constitutes the signal. 


where the detunings A c = cj 0 i — cj c , A p — — uj p and we have taken a p = iQ p , with a real Q p . We will also 

denote the coupling operators as Lij = y / T \j \i) (j\ where Tij is the decay rate from j th energy level to zth level. 
The above Hamiltonian is valid for an atom with one input-output port. Such a setup can be achieved by placing 
the transmon at the end of a semi-infinite transmission line. With the transmon initially in the ground state 10), 
the probe field does not interact with the 1 — 2 transition. The probe field is scattered only when the control field 
has a photon which excites the transmon to the state |1). By monitoring the output probe field using homodyne 
detection the presence of the single control photon can thus be inferred. 

The single photon in the control field can be modeled as the output from a fictitious cavity with damping rate 
K,(t). By modulating K,(t) = £(t)/y / f t °° |£(s)| 2 ds, the shape of the temporal wave packet £(t) containing the single 
photon can be arbitrarily set [38]. The master equation of the setup is 

p — —i [H, p] + (T> [Loi] + 'D [L12]) p + K,(t)T) [a] p — yjhz(t) C [a, L01] p 
= ^ca vP: (5) 


where p is the total density matrix of the source cavity and the transmon, a (a’*’) is the annihilation (creation) 
operator for the photons in the cavity and the dissipation super-operator is given by T> [c\ p = cpc ^ — |c^cp— |pc^c. 

We have also defined a Liouvillian C cav and a coupling super-operator C[ci,C 2 ]p = ^4,cipJ + cj,C 2 J for 
shorthand. 

The output probe field that is measured using the homodyne detector is given by the standard input-output 
formalism as = Q P + T 12 . However, we are only interested in the change of the probe field with and without 

a single photon in the control field, i.e. only on L 12 . The interesting part of the homodyne current can thus be 
defined as 

j(t)dt = y/rj (e^L^it) + L 2 i(t)^ dt + dW(t ), ( 6 ) 

where 0 < 77 < 1 is the efficiency of the homodyne detector and <j) is the phase of the local oscillator that specifies 
the quadrature of the measurement. dW{t) is a Wiener increment with a mean E[dW(t)\ — 0 and variance 
E[dW(t) 2 ] = dt where E[-] is the ensemble average. To convert the above time trace of the homodyne current 
into a binary flag, we define a signal S = f* f j(t)f(t)dt , where t m = tf —U is the measurement time window and 
f(t) is a linear filter function. We will initially take f(t) to be just a square pulse with value 1 when U < t < tf 
and 0 otherwise. 

To characterize the performance of this setup as a photon detector, we can define the signal to noise ratio (SNR) 
as 

SNR = .M , (7) 

v/Var^i] + Var[S'o] 

where ,S'q/i is the signal with 0 or 1 photon in the control field and Var[X] = E[X 2 ] — E[X] 2 is the variance. The 
SNR can be calculated from the above master equation (5). The average and variance in this case are 
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( 8 ) 


E[S 0 ] = 0 
E[S%] = t m 

E![Si] = ^/rjJ ( y)dt 
E[Sl] = dt 2 E\j(ti)j(t 2 )\. 

The two time correlation function can be evaluated using the quantum regression theorem [39] as 

E\j(ti)j(t 2 )} = 0(t 2 - ti) (^Tr ((e^Lia + e~^L 21 )T(t 2 - h )(e i *L 12 p(t 1 ) + pit^L^j) + 5(h - ti)) 

+ 0(ti - h) (pTr ((c^Lia + e~^L 21 )T(t i - t 2 )(e^L 12 p(t 2 ) + e~^p(t 2 )L 2 1 )) + $(ii - i 2 )) (9) 

where T(t 2 — ti)V(ti) = ¥(t 2 ) with ¥(t) = e l< ^Li 2 p(t) + e~ l( ^ p{t)L 2 i. The time evolution operator T{ti — £ 1 ) is 
evaluated by solving ¥ == C cav ¥. By definition, the step function 0(£) = 0 for t < 0 and 1 otherwise. 

The signal to noise ratio is an ideal measure for Gaussian statistics but is not necessarily a good measure for other 
distributions that need higher orders moments for complete description. In order to collect such statistics of the 
signal distribution and define more relevant measures, we use the formalism of stochastic master equations(SME). 
These equations describe the evolution of the system under measurements and can be considered as an unravelling 
of the average system dynamics described by the above master equations [40]. With the tunable cavity as the 
photon source, the SME is given by 


dp = £ ca vpdt + ^/rjM [A i2 ] p dW(t ), (10) 

where we have defined a measurement super-operator M [c] p — {e' l ^cp-\-e~ l( ^pc^) — (e^c + e - ^c^) p, that describes 
the back-action of the measurement on the evolution of the system. 

To get the distribution of the signal, we numerically solve the above equations with n = 0 and n = 1 photon 
in the control field. Each run of this simulation is called a trajectory which gives a particular value of S. The 
number of trajectories should be high enough to get a probability distribution for S. In this case, we can define 
fidelity of single photon detection as 

F = P(0\S < S£)P(S < ScO + P(1\S > ST)P(S > SH, (11) 


where P(0|5' < So) is the conditional probability to have 0 photons in the control field given that the measured 
integrated current S was found to be less than or equal to a predefined threshold Sq . P(l|£ > Si) is defined 
analogously. 

As first shown in [35], the signal generated using a single transmon in the above setup cannot overcome the 
quantum noise. This can be seen from the distribution of the signal shown in Fig.3, which is obtained from the 
numerical simulation of the above stochastic master equations. The distribution is obtained from 2000 trajectories 
for n — 0 and n — 1 control photons, assuming no additional losses such as dephasing in the transmon and with 
perfect homodyne detection (rj — 1). In this simulation, we have also taken the single photon to be in a Gaussian 
wave-packet given by 


m 


r^V /4 f-r 2 ph (t-T ph ) 2 

-2V ) ° XP (-4- 


where T p h gives the width of the wave packet and T p h is the time of arrival of the peak of the wave packet to the 
transmon, with the normalization condition f \^(t)\ 2 dt — 1. The parameter range is chosen to give optimal SNR 
which was calculated to be around 0.7 and correspondingly the fidelity was found to be around 70%. 

The reason for SNR < 1, can be understood heuristically as follows. A single atom can process only one 
excitation of a transition per interaction time. Even if we consider that the control photon is completely absorbed 
by the atom, it can only scatter a single photon from the coherent probe field. As the noise of the coherent field 
is of the order of a photon, this displacement in its amplitude is not enough to distinguish between the cases of 
having n — 0 and n — 1 control photons. In fact, it can be seen that driving the 1 — 2 transition with the probe 
field actually reduces the excitation probability of the 0—1 transition (i.e. the incoming single photon is not 
completely absorbed by the atom) and thus we don’t even reach the limit discussed above (See Fig.4). 
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Figure 3. Histogram of the integrated homodyne current S with(red) and without(blue) a photon in the control field, using a single 
transmon (N = 1) at the end of a transmission line. The number of trajectories in both the cases (n = 0 and n — 1 control 
photons) was 2000. The parameters used in the units of Toi are Aoi = A 12 = 0, Q p = 0.35 and V \2 = 2 while the quadrature of 
measurement is set with (f) = 7t/2. The input control photon is of Gaussian temporal shape with T p h = 0.8 and T p h = 4. The time 
of measurement t m was optimized to give the best SNR. As is evident from the distribution, the SNR in this case is less than 1. 




Figure 4. Excitation probability of the 0—1 transition of a transmon: In the absence of a coherent probe field on resonance with 
1 — 2 transition, the single photon in the control field is completely absorbed by the transmon, if we enclose it in a rising exponential 

wave-packet £(t) = ©(£ — T p h) -\J^ph exp ^ P ^ 1 ^ with F p h = Tot. Switching on the probe field, reduces the maximum excitation 
probability as can be seen in the second panel. 

3.2. Beyond a single transmon 


In order to overcome the limitations as discussed above, we look at using more than 1 transmon to displace 
the coherent probe over the limit set by the quantum noise. As a first step in this regard, we briefly review a 
convenient tool known as the (S,L : H) formalism [41,42] to derive the master equations for connected quantum 
systems. In this formalism, each subsystem is described by a triplet G = (S,L,H), where S is the scattering 
matrix, L is the vector of coupling operators and H is the Hamiltonian of the subsystem. Once the triplets are 
identified for each of the subsystem, the total triplet for the composite system can be written using the following 
products. 

The series product < of the triplets describes feeding the output from one subsystem into another 


G 2 <G 1 


^2ft,^ 2 T 1 +L 2 ,if 1 +if 2 + - 



( 12 ) 


The concatenation product S3 is used for composing subsystems into a system with stacked channels 
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Figure 5. Schematic setup of N transmons coupled to a transmission line. Each of the transmons is probed by a coherent field and 
the change in the reflected or transmitted amplitude is measured by a homodyne detector. 


G 2 BB Gi = 



,tf 2 + iA 


(13) 


The (S', L, H) formalism can also handle feedbacks. However, as we won’t use them in this particular review 
we refer the readers interested in feedbacks to the above references. 

Using the above defined products, we can write down the (S, L, H) triplet for the whole system 


(14) 



( 


\ 

Gtot = 

Stot 1 


5 Htot 


l 

\Lnj 

/ 


from which we can extract the corresponding master equation as 

n 

p= -i [H to t,p] + Lj V [Li] p- 


(15) 


The output from the i th channel is just given by L*. The (S, L, H ) triplet for the case of a single transmon with 
three levels in front of a mirror is 


Gtr — I f-2 



(16) 


with the same Hamiltonian as in the previous section. The (S, L, H) triplets for the tunable cavity and the 
coherent probe in their corresponding rotating frames are 

G cav = (1, y/ft(t)a, 0) (17) 

and 

G« p = (l,a p ,0). (18) 

With these triplets, we are now ready to look at different approaches with which we can connect different 
transmons. A first approach that is experimentally straightforward to implement is to couple many transmons 
to a transmission line one after the other (see Fig. 5). One can imagine two different scenarios in this case. 
The distance between the transmons could be much less than the wavelength of the incoming fields or it could be 
several factors higher. Both of theses scenarios were analyzed in [35] and found to be not helpful for the problem of 
single photon detection. The heuristic understanding for these is as follows. In case the transmons are closer than 
the transition wavelengths, the energy levels hybridize to form a larger atom with a different normalized coupling. 
In this case we are back to the case of a single transmon, albeit with a different set of parameters. While this 
limitation could in principle be overcome by separating the transmons over larger distances, a Kramers-Kronig 
type relation between the interaction and losses prevents any useful gain in the signal. To go beyond these 
limitations, we need a unidirectional coupling between the transmons i.e. a cascaded setup as discussed in [36]. 
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The setup considered is schematically shown in Fig. 6. We achieve unidirectional flow of the fields by placing 
the transmons at the end of transmission lines (similar to having atoms in front of mirrors) and connecting them 
using microwave circulators. The master equation for this setup is derived as follows. Following the product rules 
in Eqs. 12 and 13, the (S,L : H) triplet for the setup consisting of N cascaded transmons with a tunable cavity 
as the photon source and a coherent probe of strength a p is 


Gtot = G{? < g£ } <a (G cav ffl G ap ) 

I \/n(t)a + Aqi 


= la, 


y OL p + A12 




where 


Htat = 


N i 

RtP + — y/K,(t) ^Aion — a^Aoi^ + — (a p A2i — a*Ai2) 

3 = 1 1 1 

N ( N \ N ( N N 

E Z + ±E E (l&'L® - 

3 = 1 \k=j +1 ) 3=1 yfc=j + l y 


+ 2 i 


(k) 

and the collective operators A ij = ^ Lj • . This gives the master equation (after some algebra) as 

k =i 

N 


p= -i [H eS , p] + ^2 {p [ i oi ) ] + v [^l^]) P + K(t)V [a] p - y/n{t) C [a, A 0 i] p 

3 = 1 

-E E (c[LS»,4?]+c[ig,iS 


j=l k=j+1 


Ip 


(19) 


( 20 ) 


( 21 ) 


where the effective Hamiltonian is if e ff = ^2 with |0) (0|^ + |2) (2 |^ +Q P (L^ TT^). 

k=1 

We have once again chosen the normalization for the probe held such that a p = fl p e Z7r//2 , with Q p being a real 
number. A master equation using the Fock state formalism could also be derived similarly. 

We once again numerically simulate the above equations and corresponding stochastic master equations to 
calculate the signal-to-noise ratio and the fidelity of photon detection. The main results are shown in Fig. 7. As 
we can see from the figure, the unidirectional coupling helps accumulate the effects from each of the transmon and 
we can break-even the noise limit with N = 2 transmons. The resulting SNR can be fit to a simple y/~N curve, 
which would be the expected behavior if we considered each of the scattering events to be independent of one 
another. The results in Fig. 7 are once again for the ideal case of no dephasing in the transmons and with perfect 
homodyne detection. The setup also seems to work if we take into account some deviation from the ideal setup, 
though with reduced fidelities. We will not reproduce these results here and would refer the reader interested in 
full details to reference [36]. 


Transmon 



Figure 6. Schematic setup of cascaded transmons for QND detection of microwave photons. Placing the transmons at the end of 
transmission lines and connecting them via circulators makes this a cascaded system where the fields travel unidirectionally. 
































Figure 7. Signal to noise ratio (SNR) and the fidelity of photon detection F as a function of N , the number of cascaded transmons 
using the Fock state master equation. Both the control and the probe fields were taken to be on resonance with all of the transmons. 
The coupling of the individual transmons were tuned to optimize the SNR (Refer [36] for full parameter list). The input photon 
was taken to be in a Gaussian wave packets. The insets show the signal distribution for N = 2 and N = 8 transmons. Note that 
for the case with an input control photon n = 1, the distributions deviate from a normal distribution with the increasing number 
of transmons, suggesting a memory effect. In such cases, we believe the fidelity of photon detection becomes a better measure than 
SNR. 



Circulator 


Figure 8. Schematic setup of a transmon with a cavity for the probe field. The setup can be extended by adding more units 
(transmon + cavity) using circulators to get a cascaded setup. The field reflected from each cavity is measured using a homodyne 
detector to infer the presence of the control photon. 

One of the main concerns in experimentally implementing the above setup is the need of circulators. Currently, 
these are bulky devices that are lossy and are off-chip. On chip circulators [43-45] or other unidirectional wave¬ 
guides (a la the quantum hall edge channels [46]) will make this scheme more attractive to implement. The 
performance of these setups could be further improved by using a cavity for the probe as shown in [37]. We will 
briefly review this in the next section. 


3.3. Back to cavities 


As mentioned in the introduction, early experiments in QND detection of microwave photons were done in 
cavity QED. A cavity field can interact with the atom over several cycles, as the interaction time between the 
atom and the cavity field is much less compared to the cavity lifetime. However, if we are interested in detecting 
an incoming photon from arbitrary source (say from a different quantum node), we have to capture the photon 
in the cavity. To fully absorb the photon, the cavity needs to be of wide bandwidth and with strong coupling to 
the transmission line. High coupling however reduces the cavity lifetime, which degrades the performance of the 
photon detection schemes. However, as the probe field is only an auxiliary field for detecting the incoming control 
photon, one could imagine having a cavity for the probe only while the control photon is still an itinerant one. 
Such a setup (Fig. 8) was proposed and analyzed in [37]. 
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We will first consider a single unit comprising of a transmon with a cavity for the probe field. In the rotating 
frame of the input fields, the Hamiltonian can be written including the coherent drive as 

H = 6i\l)(l\ + (6i + 5 2 ) 1 2) (2| - iE(b - 6+) - ig(ba 21 - b'a 12 ) (22) 

where b is the annihilation operator of the probe cavity, E is the amplitude of the drive and g is the coupling 
strength between the cavity and the 1 — 2 transition. The stochastic master equation for this unit, using a cavity 
as a source of the control photon can be written as 

dp = -i [H , p\ + (V [Loi] + V [Li 2 ]) p + K a (t)V [a] p + K b V [b] p - ^K a (t) C [a, L 0 1 ] p. (23) 

The field output from the probe cavity is measured using a homodyne detector and the corresponding homodyne 
current is 

j(t)dt = Vwj (e~ i4, b + e i4, b f ) dt + dW{t). (24) 

Similiar to previous section, the performance of this setup was analyzed numerically and the fidelity F was shown 
to improve from 70% to 84%. This was achieved with the help of the cavity and by using an optimal linear filter 
f(t) that takes the form of the expected homodyne current when the control field has a photon. The results 
also assume that the decay from the state |2) of the transmon into the transmission line can be suppressed to 
T 12 = O.lToi. Lifting this restriction and using the usual transmon limit with T 12 = 2Toi, the fidelity was found 
to be 81%. Further improvements were shown to be possible by cascading the units similar to the previous section. 


4. Summary 

Microwave quantum optics using superconducting circuits, dubbed circuit QED, is an alternate approach for 
studying light-matter interaction that has been developing rapidly over the last few years. By confining the 
electromagnetic field to 1-dimension, strong coupling regime can be reached in these setups that has already led 
to plethora of interesting results. A single photon detector in microwave regime will fill a major gap that would 
make these setups even more attractive for studying quantum optics and for developing quantum information 
technologies. The proposals we discussed in this article are aimed towards this goal. Especially, a QND single 
photon detector would lead to several interesting applications as discussed earlier. The setups (or similar ones) 
discussed for QND detection are close to experimental realization and will become even more attractive with on- 
chip circulators and quantum limited amplifiers. We believe that new and exciting physics in microwave quantum 
optics will be studied using single photon detectors in the near future. 
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